Synchronization universality classes and stability of smooth coupled map lattices 
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We study two problems related to spatially extended systems: the dynamical stability and the 
universality classes of the replica synchronization transition. We use a simple model of one di- 
mensional coupled map lattices and show that chaotic behavior implies that the synchronization 
transition belongs to the multiplicative noise universality class, while stable chaos implies that the 
synchronization transition belongs to the directed percolation universality class. 
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I. INTRODUCTION 

Chaos is a widely studied phenomenon. It represents 
the generic behavior of nonUnear systems and it is also 
assumed to be the microscopic mechanism for the am- 
pUfication of uncertainties about initial conditions. The 
exponential growth of an initial small uncertainty is mea- 
sured by the maximum Lyapunov exponent (MLE). How- 
ever, in high-dimensional systems, and especially in ex- 
tended systems with short-range coupling, the origin of 
uncertainty may be different from the previous one. To 
be more specific, the systems we take into considerations 
are defined on a lattice with short range interactions. 

Just to give a simple example of the different role of 
chaos in low-dimensional with respect to extended sys- 
tems, let us consider an array of uncoupled chaotic maps. 
This is clearly a chaotic system; a localized uncertainty 
will grow in time and will remain localized. If there is 
a coupling between neighboring maps, a localized uncer- 
tainty will propagate to neighboring sites. This prop- 
agation is a linear phenomenon, and is hidden by the 
exponential amplification due to the chaotic dynamics 
of maps. If we now consider non-chaotic {i.e., stable) 
coupled maps, only the spreading mechanism may be 
present. Since the individual maps are stable, infinites- 
imal uncertainties are always absorbed, but finite ones 
may survive. One of the most striking examples are 
class-3 cellular automata (CA) that are equivalent 
to superstable maps, but may nevertheless be denoted 
"chaotic" because a localized uncertainty may spread to 
the whole system. 

While cellular automata may be seen as coupled maps 
of extreme type, that only take values zero and one, 
CA-like behavior has been observed also in continuous 
maps, albeit with discontinuities or extremely sharp tran- 
sitions 2]. In these cases the maximum Lyapunov expo- 
nents (if defined) is negative but defects and information 
may propagate and amplify in a non-exponential way. 
When finite, these systems are generally periodic, but 
the period grows at least exponentially with the size, so 



that the transient is the only observable state Q, |M ^ ■ 
This behavior is denoted as stable chaos. 

As mentioned before, evaluation of the MLE is the 
most direct method to establish the stability properties 
of a system but is somewhat limited for extended sys- 
tems as coupled map lattices (CMLs). An alternative 
method for the measurement of uncertainty is based on 
the synchronization of replicas of a system |^ . Typ- 
ically, one fraction of each replica is added to the other 
and the distance between them vanishes as this fraction 
varies 0, 1^ • Another mechanism involves the addition 
of the same external noise to both replicas 9]. In both 
cases, there is a control parameter p that measures the 
strength of the synchronization. 

We may have different scenarios, according with the 
degree of unpredictability of the system. Chaotic systems 
are expected to amplify the distance between replicas. 
For a value oip slightly below the synchronization thresh- 
old, some patches may synchronize for some time, after 
which they will separate. This picture resembles that of 
a growing interface that may stay pinned to local traps. 
From field theory studies, such a behavior is denoted mul- 
tiplicative noise (MN) and is equivalent to the behavior 
of the (bounded) Kardar-Parisi-Zhang equation, which 
describes the behavior of a growing surface that tends 
to pin and is pushed from below [li IH IH El El. 
On the other hand, stable systems have a negative MLE. 
So, replicas should naturally synchronize once their dis- 
tance is (locally) below the threshold of validity of linear 
analysis. However, when the (local) difference is large, 
non-linear terms may maintain or amplify this distance. 
In this case synchronized patches may be destabilized 
only at the boundaries. Again, theoretical studies as- 
sociate such a behavior to that of directed percolation 
(DP) 0, 0,^^. In other words, one may illustrate this 
behavior by saying that the synchronized state is simi- 
lar to the void state of the field theory descriptions, and 
that this void state is absorbing in the second case and 
unstable in the first case |l3 . 

It has been noted that CMLs with continuous chaotic 
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maps lead to synchronization with a MN character while 
the synchronization of CMLs where the local map is dis- 
continuous 0, ^3 , or with very sharp transitions Q , be- 
long to the DP universality class. The use of simple mod- 
els where the transition between the two behaviors may 
be triggered by a variation of control parameter has been 
addressed for the stochastic synchronization of CMLs in 
Ref. ^3 and in growth models in Refs. [l3Lll3 |. Random- 
ness can be safely used if the system is chaotic although 
the systems under study are deterministic. This is not 
the case when there is stable chaos, unless coarse graining 
is used. 

In what follows we investigate the relation between sta- 
bility and the synchronization character of CMLs. We 
present a model that, by changing continuously a single 
parameter, can vary smoothly from chaotic to chaotically 
stable (non chaotic but unpredictable) behavior, and fi- 
nally to cellular automata behavior. From the point of 
view of the synchronization transition, the model exhibits 
a transition from the MN universality class to that of DP. 
The transition from chaos to stable chaos occurs at the 
same value of the parameter as the transition from the 
MN to the DP universality class. Moreover, we show that 
these transitions occur for smooth and continuous maps, 
and that at the transition stable synchronized trajecto- 
ries are still present. The transition from standard chaos 
to stable chaos occurs through the appearance of tran- 
sient chaos, i.e., the system in the stable chaotic phase ex- 
hibits a transient chaotic behavior before being attracted 
into the stable-chaos trajectory. 

The paper is organized as follows. In Sec.^we present 
the CML model and the synchronization mechanism. In 
the following section we present the main numerical re- 
sults and in Sec. we end with some conclusions and 
directions for further work. 



II. THE MODEL 



A coupled map lattice F defines a fiow on [0,1] , 



„t+i 



F{x^) with X = (xo, . . . ,xn-i) and i = 0, 1, 



It is defined locally by 

x*+i=/((l-2eH + e(x*_i+x*+i)) (1) 

with i = Q, . . . ,N — 1 and periodic boundary conditions, 
< e < 1/2 and / : [0,1] [0,1]. In what follows 
e = 1/3 so that 



(2) 



with (Ti = {xi-i -I- a;i -I- a;i+i)/3. 

The model is defined by the choice of / that depends 
on a parameter a such that for large values of this pa- 
rameter, the CML behaves as an elementary cellular au- 
tomaton (super-stable system), and in the other limit the 
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FIG. 1: (Color Online) The graph of /(s;) for three values of 
a. 

individual map is chaotic. One way to define / is by 



f{x:a) 



'(6a;)°/2 

1- [6(1/3 -x)|V2 
|6(x- 2/3)172 
1 - (6(1 - a;))V2 



< a; < 1/6, 
1/3 < a; < 1/2, 
1/2 < a; < 5/6, 
5/6 < a; < 1, 



(3) 



where 1 < a < oo. This function was introduced in 
Ref. 19] in the context of stochastic synchronization. 

The graph of / is shown in Fig. ^ and a typical space 
time pattern is shown in Fig. [3 For a choice of the 
parameter a > 1.8, the CML is asymptotically equivalent 
to a cellular automaton (CA) with triangles of zeroes and 
ones. The function / was chosen in such a way that the 
CML in the limit of large a will behave as the "chaotic" 
elementary CA rule 150 JJ. 

Chaotic systems are characterized by the divergence 
of initially near trajectories, and a quantitative indica- 
tor for chaotic motion is (generally) a positive maximum 
Lyapunov exponent (MLE) A(a;°), that depends on the 
trajectory. In order to probe the configuration space, we 
numerically evaluate Xm.Ti the average over M initial 
randomly chosen configurations a;" of finite time Lya- 
punov exponents during T time steps. 

We study the synchronization of CMLs by considering 
the following scheme: 



4^^ = H^' 
= (1 



(l-p)/(T*)+p/(a*), 

P)\fi^l)'fir!)\, 



(4) 
(5) 
(6) 
(7) 



where = (?/i_i + yi + j/i+i)/3. Synchronization occurs 



when • 



Vz, which is equivalent to a vanishing norm 
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FIG. 2: (Color Online) Space time pattern of the CML of 
Eq. ^ with a — 1.9 and A'^ = 256 drawn horizontally for a 
total of T = 300 time steps drawn vertically from top to bot- 
tom. The initial configuration x" is chosen randomly. The 
color code assigns white (black) whenever a;' = (1) and a 
rainbow color scale for other values of a;' starting with red 
for values near zero. Patches of CA behavior (rule 150) ap- 
pear after a short transient and will eventually fill the whole 
pattern. 

M* of It* defined by 

i 

The synchronization threshold is called Ps- In numerical 
simulations, we usually evaluate (u), defined as an aver- 
age over T time steps after a relaxation time trei which 
is then averaged over M random initial conditions. In 
Fig. Elwe show u* for several values of a and p^. 

The linear stability of u is measured by the transverse 
Lyapunov exponcnt(TLE) Aj^ defined by |23| 

Ai(a;") = log(l-p) + A(a;0). (9) 

The synchronized state becomes stable with respect to 
infinitesimal perturbation when A^ = so that the value 
of p at synchronization, pi, is given in the linear approx- 
imation by 

= l-exp(-A). (10) 

This does not mean that synchronization actually takes 
place for this value of p. It may happen that the system 
stays unsynchronized for p > pi. In Ref ,7,] this behavior 



has been associated to discontinuities of the single map 
taking as reference the Bernoulli shift. We shall show in 
the following that even smooth, non-chaotic maps may 
correspond to ps > pe- 



III. RESULTS 

We first discuss the stability properties of the CML 
presented in the previous section and then the character 
of the synchronization transition. If 1 < a < Oc the 
CML is chaotic and if a > Oc the CML exhibits stable 
chaos. We found Uc — 1.8142(2) as discussed below. For 
I < a < ttc almost any trajectory starting at any initial 
state a;° with a;° chosen randomly between zero and one 
has a positive MLE. However, if the initial state is such 
that = 0, 1 randomly, the CML will evolve as the 
elementary CA rule 150 with a negative MLE (actually 
A = — oo). 

The configuration space of the system is the uni- 
tary hypercubc [0, 1]^, whose vertices correspond to au- 
tomata configurations (zeroes and ones). Starting from 
one vertex, the time evolution will visit only vertices. 
There is an attractor near the vertices of the hypercube 
that brings any trajectory with initial state in this at- 
tractor toward these trajectories CA trajectories. The 
transition from chaos to stable-chaos corresponds to the 
transition where this attractor fills the hypercube. 

For a > flc the system exhibits transient chaos: ini- 
tially nearby trajectories diverge and are finally attracted 
to a CA type evolution. Numerically one observes a pos- 
itive finite-time MLE that (suddenly) becomes negative. 
This fact is used to find the value of Oc as shown in Fig.O 
For different values of a, the MLE of M orbits starting 
from randomly chosen initial states is evaluated and the 
fraction of chaotic orbits is shown as a function of time. 
The first value of a for which this fraction is not one is 
taken as flc. For a > Gc the fraction of chaotic orbits 
goes to zero as i ^ oo so the measure of the set of initial 
states showing stable chaos goes to one. 

Let us now turn to synchronization. In Fig. 0] we show 
the dependence of the order parameter (u) with p. For 
values of a close to the critical value there is a wide spread 
in the values of u as shown in Fig.0](b) making it difficult 
to define Ps- By letting the system to relax longer, we 
find a well defined synchronization threshold as shown 
in the inset of the figure. In Fig. O we show the phase 
diagram of ps vs. a. Note the large jump of ps at Oc. 

For a < Qc, Ps — 0.5 and for values of p close to this 
value, the two replicas are forced to stay a small distance 
apart, so that only the linear (chaotic) separation is ac- 
tive. The synchronization transition thus occurs in agree- 
ment with Eq. H10|l , the difference field may be contract- 
ing or expanding, according with p, but the synchronized 
state is not absorbing (panel (a) in Fig.|HJ). By increas- 
ing a, the Lyapunov exponent drops to — oo (panel (b) 
in Fig. 13 and we cross the phase boundary. By lowering 
p, the non-linear mechanism for the diffusion of the dif- 
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FIG. 3; (Color online) Fraction of chaotic orbits as a function FIG. 5: Synchronization threshold ps vs. a. Labels refer to 
of time for different values of a for a total of M = 100 random patterns m Fig. |Bl 
initial states, N = 1,024 and T = WON. 
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FIG. 4: (Color online) Dependence of the order parameter 
{u) on p near the critical value ac- In (a) T = trei = 5 x 10*, 
M = 100, and = 1,024. In (b) the standard deviation is 
rather large for the same values of T, t^ei, M and A'^ making it 
difficult to assign a value to Ps- However, by making trei ~ 10^ 
there is a well defined transition at ps — 0.303 as shown in 
the insert. 



ference field is activated and sustains the differences be- 
tween the replicas. The cluster of non-synchronized sites 
is now connected, like in DP clusters (panel (c) in Fig.O. 
By further increasing a, the linear mechanism becomes 
quickly contracting, and the non-synchronized cluster is 
more reminiscent of cellular automata ones (panel (d) in 
Fig. El. 

By extensive numerical simulations we found the values 
of the critical exponents S, [3 and z defined by 



with g an unknown function. For a below the criti- 
cal value ttc the synchronization transition belongs to 



the bounded Kardar-Parisi-Zhang (BKPZ) universal- 
ity class ( S = 1.10 ± 0.12, P = 1.50 ± 0.15, and 
z = 1.53 ± 0.07) and for a > Oc it belongs to the directed 
percolation (DP) universality class {S = 0.159464(6), 
(3 = 0.276486(6) and z = 1.580745(1) as shown in Fig.|7| 
The precision of the reference values of the critical ex- 
ponents refiect the difficulties in measuring the critical 
behavior of systems in the presence (absence) of absorb- 
ing states. Numerical experiments concerning DP are 
much more precise that those concerning MN, due to the 
absorbing character of the void state. In both cases the 
most sensitive measure concerns the spreading rate of 
a defect, but while in some DP simulations (not in this 
case) it is sufficient to simulate the behavior of the system 
around the spreading damage (the evolution of the void 
state is trivial), in MN experiments one has to simulate 
the whole system. This implies that these DP experi- 
ments immediately reach the asymptotic limit for what 
concerns the spatial dimension, while MN experiments 
suffer from finite-size spatial effects. Moreover, DP sys- 
tems may be discrete, while MN ones are simulated using 
floating-point arithmetic. 



IV. DISCUSSION 

Our results suggest that the jump from one universal- 
ity class to the other occurs at the value Uc where the 
jump from chaos to stable chaos occurs. In other words, 
the microscopic stability properties of the CML are re- 
sponsible for its macroscopic synchronization transition. 

Near the synchronization threshold, the difference field 
behavior is dominated by the spreading properties of non- 
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(a) (b) 




FIG. 6: (Color Online) Synchronization space time patterns 
near pa for different values of a. See Fig. |3 for reference, 
(a) a = 1.8 < ac (A > 0), p = 0.48 < ^^1.8): the differ- 
ence field is expanding but low in value, and the synchro- 
nized state is not locally absorbing (b) a — 1.9 > Uc (A < 0), 
p — 0.48 3> Ps(1.9): the difference field is contracting, and 
the synchronized state is absorbing, even though occasionally 
patches of non-synchronized sites appear, (c) a = 1.9 > ac, 
p = 0.36 > Ps(1.9): the difference field is sustained by the 
nonlinear term, the unsynchronized sites form a connected 
cluster, (d) a — 3.5 3> flc, p = 0.36 > Ps(3.5): the lin- 
ear contribution to the difference field is (almost) absent, the 
non-synchronized cluster is connected and has more "holes" 
than in case (c). N ^ 256, T = 256. 



synchronized patches, which is exactly what is generaUy 
measured for the determination of the critical properties 
of the transition. The evolution of the difference field, 
which is in average small at the transition, is related 
to the linear stability of the individual map: systems 
which are linearly expanding (positive MLE) exhibit a 
difference field which is locally small everywhere, with 
occasional fluctuations. In this case one has a local com- 
petition between the synchronization strength and the 
diverging character of the dynamics. Systems exhibiting 
stable chaos, however, are quite different. The difference 
field cannot be small, otherwise it would spontaneously 
vanish due to the negativity of MLE. Therefore, it has 
to be large in some patches of small extension (since at 
the transition it is small in average), and is maintained 
and propagated by the nonlinear part of the dynamics 
that amplifies "instantaneously" the differences between 
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FIG. 7: (Color online) The values of the critical exponents S, 
13, and z for different values of a. The vertical line is drawn 
at ac = 1.8142. 



replicas. This amplification means that even if two state 
values in a certain location of the two replicas are simi- 
lar, they may become very different in just one step. This 
effect can obviously be generated by the explicit discon- 
tinuities in the local evolution rule [3| , but may manifest 
also in continuous maps, as shown by our toy model. 
Actually, our model exhibits this behavior when the dy- 
namics has reached the "cellular automaton phase", i.e., 
after the eventual transient chaotic regime. In this CA 
phase the effective dynamics is discontinuous, but this 
behavior is not easily recognizable just by looking at the 
local evolution rule. Notice that DP-like transition can 
be exhibited also by truly chaotic systems with disconti- 
nuities, like Bernoulli shifts 0. Our results are consistent 
with recent findings using a stochastic field theory that 
recovers both MN and DP behavior in a unified way re- 
lated to the change of stability of the system p^ . 



V. CONCLUSIONS 

We have investigated the relation between the type of 
unpredictability and synchronization transition univer- 
sality classes, using a coupled map lattice model where 
the local dynamics is given by a continuous map that 
can exhibit chaotic or stable character. This model may 
exhibit either chaotic unpredictability or stable chaos. 

We have analyzed the synchronization properties of 
two replicas, coupled asymmetrically with a parameter 
p. Linear stability analysis gives an estimate synchro- 
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nization threshold p(, related to A. This is actually the 
case for smooth chaotic maps. 

We have shown that this transverse direction may stay 
expanding (i.e., the system does not synchronize) for 
p > p£. In this case the synchronized trajectories are 
actually absorbing, but the measure of initial conditions 
that bring toward them is vanishing. 

We expect to observe similar behavior in other, more 
physical systems, and in particular in biological ones: 
the kind of unpredictability (chaos) observed in biologi- 
cal systems cannot be ascribed to a chaotic behavior of 
cells, that in viable conditions are rather stable unities, 
but is presumably due to the coupling among them, in a 
way that is reminescent of cellular automata "chaos" . 

Actually, cellular automata have become a widely 
adopted modeling tool, due to the simplicity of descrip- 
tion of the interaction rules and the efficient simulat- 
ing techniques. They have been used to model a wide 
variety of systems, from physics to engineering to biol- 



ogy. However, one may ask where is the transition from 
the continuous, microscopic world to the discrete, macro- 
scopic description of cellular automata. This is exactly 
the question addressed in this paper, and we have shown 
an example in which this transition may be located quite 
precisely. 
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